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Magnetic relaxation and dipole-coupling-induced magnetization 
in nanostructured thin films during growth: A cluster Monte Carlo study 

R. Brinzanik0 P.J. Jensen, and K.H. Bennemann 

Institut fur Theoretische Physik, Freie Universitat Berlin, Armmallee 14, D- 14-195 Berlin, Germany 

(Dated: February 1, 2008) 

For growing inhomogeneous thin films with an island nanostructure similar as observed in experi- 
ment, we determine the nonequilibrium and equilibrium remanent magnetization. The single-island 
magnetic anisotropy, the dipole coupling, and the exchange interaction between magnetic islands 
are taken into account within a micromagnetic model. A cluster Monte Carlo method is developed 
which includes coherent magnetization changes of connected islands. This causes a fast relaxation 
towards equilibrium for irregularly connected systems. We analyse the transition from dipole cou- 
pled islands at low coverages to a strongly connected ferromagnetic film at high coverages during 
film growth. For coverages below the percolation threshold, the dipole interaction induces a col- 
lective magnetic order with ordering temperatures of 1 - 10 K for the assumed model parameters. 
Anisotropy causes blocking temperatures of 10 - 100 K and thus pronounced nonequilibrium effects. 
The dipole coupling leads to a somewhat slower magnetic relaxation. 

PACS numbers: 05.10.Ln, 75.10.Hk, 75.70.Ak, 75.75.+a 



I. INTRODUCTION 



The investigation of low-dimensional magnetic nanos- 
tructures has become a very active field of current 
research^ The controlled preparation of different nanos- 
tructures allows for the investigation of a variety of inter- 
esting magnetic propertiesiSiiii The dependence of these 
properties on the nanostructure and magnetic interac- 
tions in such systems is still not well understood. So far, 
no consistent theoretical analysis has been performed for 
the magnetic behavior of an ultrathin film during growth, 
ranging from an island-type structure to a smooth film, in 
particular the influence of structural disorder on the mag- 
netic properties. In this study, we present Monte Carlo 
(MC) calculations of the nonequilibrium and equilibrium 
magnetization of growing inhomogeneous films. Differ- 
ent coverages below and above the percolation thresh- 
old and different magnetic interactions are taken into ac- 
count. This analysis needs a newly developed cluster MC 
method including coherent rotations of neighboring mag- 
netic islands. 

To illustrate the problem, we discuss for instance an 
ultrathin Co film grown on a Cu(OOl) substrate££ This 
system shows a perfect layer-by-layer growth for cov- 
erages 6 > 2 monolayers (ML). Here, the exchange 
coupling results in a large Curie temperature Tq > 
350 K. Below the measured percolation coverage of Op = 
1.7 ML, randomly positioned Co islands with a consid- 
erable admixture of Cu atoms are observed, exhibiting a 
comparably strong remanent magnetization up to tem- 
peratures as large as ~ 150 K. At 9p, a jump of Xc 
of about 100 K occurs^ The important question arises 
which mechanism causes the strong remanence of an 
island- type film for coverages © < ©p. On the one hand, 
this could be induced by nonequilibrium blocking effects 
due to single-island anisotropies which impede magnetic 



relaxation. On the other hand, an equilibrium magne- 
tization can originate from long-range magnetic interac- 
tions. We like to discuss these two limiting cases. 

An ensemble of magnetically isolated single-domain is- 
lands behaves like a supcrparamagnct. The single-island 
anisotropy causes a finite time-dependent magnetization 
below the nonequilibrium blocking temperature 71,. For 
thin Co/Cu(001) films this temperature is estimated to 
be Xb w 5 K^. 1 following from the Arrhenius-Neel ansatz 

?b = NK/ ln(r m T )2& Note that an island size dis- 
persion will influence the relaxational behavior and 7b 
sensitively. Second, a finite magnetization may also orig- 
inate from a collectively ordered state in thermal equilib- 
rium. Such a magnetic state, not necessarily a collincar 
one, of an ensemble of isolated islands results from long- 
range magnetic interactions^ We consider here the mag- 
netic dipole coupling between islands^ The correspond- 
ing ordering temperature Tq should be comparable to the 
average dipole energy per island which is, however, diffi- 
cult to determine for an irregular system. Assuming two 
neighboring disk-shaped Co islands containing N atoms 
each, the dipole energy E^ip of this island pair is propor- 
tional to VN. For N = 1000, we estimate E dip « 6 Ki34 

Note that both simple estimates for the tempera- 
tures 7b and 7q disagree with experimental observations. 
Hence, improved calculations are needed taking better 
into account the inhomogeneous film structure charac- 
terized by varying island sizes, shapes, and positions. In 
particular for coverages close to Op, the island coagula- 
tion leads to a larger effective island size N e g and thus 
to a larger average dipole energy E^ p cx y/N c s- Corre- 
spondingly, also the blocking temperature 7b °c N e gK 
due to anisotropy will be larger than the single-island 
estimate 7b ~ 5 K. Thus, for a disordered film struc- 
ture one may expect a remanent magnetization at much 
higher temperatures than obtained by these simple esti- 
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mates. 

For the investigation of the magnetic relaxation, the 
magnetic anistropy and the dipole interaction are taken 
into account. In case of coagulated islands, the ex- 
change coupling between islands has to be considered 
as well. For three-dimensional systems of interacting 
magnetic particles, either longei™ as well as shorter- 
relaxation times as compared to noninteracting ensem- 
bles have been calculated. Experimentally, longer relax- 
ation times for increasing interparticle interactions have 
been measuredii2i±£ The existence of a collective, spin- 
glass-like magnetic ordering was discussed^ For two- 
dimensional systems only few investigations have been 
reported, also indicating longer relaxation times for in- 
creasing interaction strengths 

In this paper, we report on MC calculations of the 
nonequilibrium and equilibrium remanent magnetization 
for nanostructured thin films during growth as function 
of coverage, temperature, and MC time— ai 7 ^ By ap- 
plication of a modified Ising model, which allows one to 
take into account magnetization dynamics, the blocking 
as well as the ordering temperatures are determined. Of 
particular concern is the consideration of structurally in- 
homogeneous systems ranging from isolated islands to 
smooth ferromagnetic films. Numerical simulations are 
unavoidable since the low symmetry of these systems and 
the complicated nature of the involved magnetic interac- 
tions preclude analytical approaches. However, the ap- 
plication of the common MC technique runs into a se- 
vere problem. Assigning a "super" spin to every island 
magnetic moment (Stoner-Wohlfarth modeli&), conven- 
tional single-spin-flip algorithms yield an extremely slow 
and unrealistic relaxation towards equilibrium for cover- 
ages where the islands are partly coagulated— Hence, 
we apply a cluster- spin- flip algorithm, which includes si- 
multaneous rotations of magnetic moments of connected 
islands— i2i We emphasize that with this method the re- 
laxation behavior and the equilibrium magnetization are 
calculated efficiently. For all film coverages, the cluster- 
spin-flip method enables an appropriate analysis of the 
influence of the anisotropy and dipole coupling besides 
the dominating exchange interaction. 

The film growth, the micromagnetic model for the cal- 
culation of the magnetic properties, and the cluster MC 
method are described in Sec. [H] In Sec. II I II we test the 
cluster algorithm for the obtained island nanostructure 
and present results for the remanent magnetization, as 
well as for the blocking and ordering temperatures. A 
conclusion is given in Sec. IIVI 



II. MODEL AND MC SIMULATION 

A. Growth mode and micromagnetic model 

For the simulation of the island-type growing film, we 
use the simple solid-on-solid Eden model— ^ Within 
this model, each additional atom is deposited on an 



island perimeter site i with probability p(Qi,Zi) oc 
cxp[A(zi) yfqi], where qi is the local coordination number 
and z% the layer index. A bilayer island growth mode is 
assumed yielding an island structure similar as observed 
for epitaxial Co/Cu(001)£ A 2 x 500 x 500) fcc-(OOl) 
unit cell with lateral periodic boundary conditions is ap- 
plied. The island density is p = 0.0025 islands per site, 
resulting in Z = 625 randomly distributed islands in the 
unit cell. For the ratio of binding parameters^ 7 - we use 
A(l)/A(2) = 0.989. 

For the obtained atomic structures, a micromagnetic 
model for the total (free) energy of a system of interacting 
magnetic islands is applied: 16,17 

E = - l -Y,L l A®)l*j{®,T)S l S J 

i>j 

+ E ^'y^ VI SiSj 3 (r ijSi ) (r ijSj )} 

i>j v 

-Y / N l (Q)K t (e,T)(Sf) 2 , (1) 

i 

with © being the film coverage and T the temperature. 
Each magnetic island with Ni(Q) atoms is treated as a 
Stoner-Wohlfarth particle 19 with a single giant magnetic 
moment /ii = /x a t rrii Ni , whose direction is confined to 
the film plane, where /z a t is the atomic magnetic moment. 
The unit vector Si — Hi/pi characterizes the magnetic 
moment direction of the ith island. The first term in 
Eq. represents the magnetic domain wall energy be- 
tween connected islands, with Ly the number of bonds 
between islands i and j, and 7-y the domain wall energy 
per atomic bond. The second term is the long-range mag- 
netic dipole interaction between the island magnetic mo- 
ments fii where = — Tj \ is the distance between the 
centers of islands i and j. The point-dipole energy is cal- 
culated by applying the Ewald summation technique over 
all periodically arranged unit cells of the thin film— ^ 
The last term denotes the uniaxial in-plane anisotropy 
energy with K% the anisotropy per atomic spin. Due to 
this anisotropy we allow for only two stable directions for 
each island moment (Sf = ±1). Thus, our system refers 
to a modified Ising model, for which during magnetization 
reversal a possible anisotropy energy barrier is taken into 
account, hence allowing the consideration of magnetiza- 
tion dynamics. Finally, due to the finite exchange cou- 
pling J between neighboring atomic spins, the internal is- 
land magnetization mj(0, T) is taken into account within 
a mean-field approximation. This leads to temperature- 
dependent effective anisotropy coefficients Ki(Q,T) and 
domain wall energy densities 7ij(0, T), as described in 
greater detail in Ref. 17. 

Equation describes a system of dipole-coupled sin- 
gle islands at low film coverages 9 <C Op as well as a 
connected ferromagnetic film at high coverages O ^ Op. 
We point out that the transition between these extremal 
cases during the film growth is described within the same 
model. The assumption of individual magnetic islands 
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with varying interactions is a good approximation as long 
as the system is laterally nanostructured, whereas for 
smooth films (here O « 2.0 ML) it represents an unphys- 
ical discretization of the system. 



B. Cluster MC method 

The magnetic equilibrium and nonequilibrium prop- 
erties are calculated by performing kinetic cluster MC 
simulations. Especially, close to the percolation cover- 
age Op, most of the magnetic islands are connected to 
neighboring islands and form large but still finite clus- 
ters. A single-spin-flip (SSF) algorithm 22 for such an ir- 
regular atomic structure yields a very slow relaxation to- 
wards thermodynamic equilibrium, since subsequent flips 
of island magnetic moments in this cluster of connected 
islands, as considered by SSF updates, are strongly hin- 
dered by the exchange energy, see Ref. 34. Thus, a ro- 
tation of the entire island cluster is very unlikely, and 
its dependence on dipole interaction and anisotropy is 
strongly underestimated. For an improved simulation of 
the magnetic relaxation, a coherent or simultaneous ro- 
tation of the spins in these clusters has to be taken into 
account. For this purpose, we propose a cluster-spin-flip 
(CSF) algorith m 2 ^ 21 in the present study. 

In a first step of each MC update, a cluster C v consist- 
ing of v connected islands is constructed by the following 
scheme: 

(a) Choose randomly a single island i, representing the 
first (smallest) island cluster C\ — {i}. 

(Jo) Add a random second island j which is connected 
to island i (Lij ^ 0), forming the second island cluster 
C 2 = {i,j}- 

(c) Construct subsequently larger island clusters C v by 
adding a randomly chosen island to the preceding cluster 
C„_i, provided that this island is connected to at least 
one of the v — 1 islands of C v -\. 

(d) Continue this construction procedure till either no 
additional adjacent islands are present or if a maximum 
allowed number A max of islands in the cluster is reached. 
From this procedure, we obtain a set of A < A max island 
clusters {C\, . . . ,C\}. 

(e) Out of this set choose one cluster C v with weight 
w v for probing ££, v uj v = 1). 

A Monte Carlo step (MCS) is defined by the usual 
condition that Z islands in the system are probed. Em- 
ploying a cluster C„ containing v islands considers the 
portion v/Z of the system in a single update. To en- 
sure that probing large clusters does not dominate the 
relaxation process, we assign the weight uj v = \jv for 
choosing C v out of the set {C\, . . . ,C\}. This definition 
implies that within a single MCS no additional relaxation 
channels are opened by the consideration of island cluster 
flips ^ We emphasize that not only is the largest possi- 
ble island cluster C\ probed for flipping, but all island 
clusters out of the corresponding set are considered. The 
island moments within an island cluster need not to be 



parallel. 

In the second step of each update, all v island spins of 
the chosen cluster C v are probed for a coherent flip. The 
corresponding flip rate Y v is calculated in the usual way 
as if these v connected islands form a single large island »i£ 
From Eq. Q), the magnetic energy of this island cluster 
as function of the in-plane angle <fi is given by 

e»(4>) = ^4^- = -2^ cos 0- cos 2 , (2) 
with the reduced magnetic field 

, _ Em S%Sf [L kl jm + y, m (Axl - 2yl) /r|] 

and the total anisotropy energy of the cluster K, v — 
Efc=i NkKk- Here, the k sum runs over all spins inside, 
and the I sum over all spins outside the island cluster 
C v . We have neglected the dipole sums J2ki( x M Zto)/ r iu 
which are usually smaller than the sums Yki( x ki) 2 l r u 

The Ising-like states S'f — ±1 of C v represent either 
energy minima which are separated by an anisotropy en- 
ergy barrier, or refer to an energy maximum and mini- 
mum. The respective energy barriers for the forward and 
backward transitions are given by 

A£« = [h v ± l) 2 Ku for \hu\ < 1 , (4) 
AEl 2) =±Ah v K v for \h v \ > 1 . (5) 

The flip rate r^, of the island spin cluster C v to overcome 
AE^ is calculated from the common Arrhenius-Neel 
ansatz£& We use a constant prefactor T = 10 9 sec -1 
which determines the time unit of the magnetic relax- 
ation in kinetic MC simulations*^ The latter case is 
treated with the usual Metropolis-type rate, using the 
same prefactor^ 

The growing thin film is characterized by a large 
amount of nonequivalent lattice sites, corresponding to 
a large number of different interaction parameters. Since 
little is known about these values, we use in our sim- 
ulation averaged quantities for the magnetic parame- 
ters which are fixed as follows, using as an example the 
Co/Cu(001) thin-film system. The atomic magnetic mo- 
ments are set to /i a t — 2.0 /iB^ The domain wall energy 
7 is adjusted to give the observed Curie temperature of 
the ferromagnetic long-range order of Tq — 355 K of a 
2-ML Co/Cu(001) film£ yielding 7 = 5.6 meV/bond. 
The exchange interaction for the calculation of the in- 
ternal island magnetic ordering is set equal to J = 
7.0 meV/bond»ii For the uniaxial anisotropy, two differ- 
ent values K = 0.1 and 0.01 meV/atom are investigated. 

In this study, we determine the remanent magnetiza- 
tion M rem (6, T, t) - Yh N i S i(t) TO *( . T )l Ni of the 
growing thin film, where t is the MC time in units of 
MCS. The simulation starts from a completely aligned 
island spin state. The choice of this initial state refers to 
experiments which saturate the magnetic system by an 
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external magnetic field and determine the remanent mag- 
netization after removal of the fieldi^ We have no evi- 
dence that magnetic arrangements end up in metastable 
states during relaxation when starting from a saturated 
state. In addition, we calculate the equilibrium mag- 
netization M e q(0,T), which is obtained by averaging 
M lem (Q,T,t) over a range of 500 MCS after the system 
has become equilibrated. M rem (0,T, t) and M eq (6,T) 
are averaged over at least 20 different structural runs. 
The magnetizations are given in units of a saturated 
monolayer (i.e., = 1 ML) at T = 0. Since the finite- 
sized unit cell undergoes eventually total magnetic re- 
versals during MC probing, accidental cancellation of a 
finite M Iem and M eq during structural and temporal av- 
eraging may occur. To avoid this, we use in this study 
merely the absolute values \M Tem \ and |M eq |. 



III. RESULTS AND DISCUSSION 
A. Island growth and CSF algorithm 

Snapshots of the atomic structure during thin film 
growth, resulting from our growth model, are shown in 
Ref. 17. The resulting static atomic structure is sim- 
ilar to the one observed for the Co/Cu(001) system^ 
In the initial stages of growth, randomly located islands 
with almost rectangular shapes are obtained. With in- 
creasing film coverages, the single islands start to coag- 
ulate and form island clusters with a still finite size. By 
analysing the percolation probability using the Hoshcn- 
Kopelman algorithm, we yield a percolation coverage of 
about 0p w 0.9 MLi 25 i 26 Continued film growth leads to 
a connected thin film. In this coverage range, the system 
still exhibits a distinct irregular nanostructure. Isolated 
island clusters vanish rapidly upon further adatom de- 
position. The coverage = 2.0 ML corresponds to a 
smooth magnetic film with two closed layers. 

At first, we investigate the effect of the cluster-spin- 
flip MC method on the simulation of the remanent mag- 
netization |M rom (0, T, t)\. Here we consider only the 
exchange interaction. For a strongly connected film 
(0 P < = 1.8 ML we test whether by use of CSF 
\M lem (Q, T, t)\ relaxes into the correct equilibrium value 
|M cq (0,T)|, starting from a fully aligned state. As can 
be seen from Fig. ^a), this condition is fulfilled, since 
different maximum allowed numbers A max of islands in 
the cluster lead to the same equilibrium value as for the 
single-spin-flip MC method. The larger the number A max 
is, the slower is the relaxation. This property is caused 
by the fact that in this coverage and temperature range 
the magnetic relaxation is mainly provided by flips of 
single islands or small island clusters. As mentioned in 
Sec. Ill Bl no additional relaxation channels are opened by 
use of CSF, hence the number of single-spin-flip attempts 
becomes reduced in favor of unprobable cluster-spin-flip 
ones. Examples for the error bars are also given. The sta- 
tistical error, which is similar for all forthcoming figures, 



results mainly from averaging over different structural re- 
alizations of the unit cell and could be reduced by using 
larger unit cells. 

We point out that the main improvement of the CSF 
with respect to the SSF method is obtained for cover- 
ages < 0p, characterized by a considerable amount 
of island cluster formation, and which is very difficult to 
be studied analytically. In Fig.^b), the magnetic relax- 
ation for = 0.8 ML is depicted using different A max . 
The equilibrium magnetization M® (Q,T) should vanish 
for < 0p, since long-range magnetic interactions are 
neglected here. Due to the use of the absolute value a fi- 
nite but small |M e ° q (0, T)| is obtained in our calculations. 
The SSF algorithm exhibits an extremely slow magnetic 
relaxation toward |M c ° q |. Even after 10 6 MCS the rema- 
nent magnetization is relaxed only to |M r ° em | = 0.60. The 
reason is that this method considers very unfavorable in- 
termediate states. Already the allowance of a few coher- 
ently flipping island spins results in a much faster relax- 
ation. The relaxational behavior converges rapidly with 
increasing A max . Using the CSF with A max = 50 or larger 
the equilibrium is reached already after ~ 100 MCS. To 
obtain a fast equilibration, the closer the coverage to 0p 
the larger the value chosen for A max . The CSF algo- 
rithm leads to a much faster equilibration also for cov- 
erages > 0p. In this coverage range, island clusters 
are still present which have only weak links to other clus- 
ters. In the following investigations, we set A max equal 
to the number of single islands Z — 625, except for cov- 
erages 3> 0p where we yield a better performance for 
Amax =100. 

Hence, by use of the SSF algorithm the exchange inter- 
action grossly dominates the MC simulations for strongly 
inhomogeneous systems. This is avoided by applying 
CSF, allowing thus for the investigation of the effect of 
the much weaker anisotropy and dipole interaction. 



B. Effect of interactions 

First, we study the combined effect of the dipole and 
the exchange interaction on the film magnetization for 
coverages < 0p. We determine equilibrium proper- 
ties which within our model are not influenced by the 
anisotropy. In Fig. |2Ja), we present \M iem (Q,T,t)\ for 
different temperatures T as a function of MC time t. 
The coverage is assumed to be = 0.8 ML. Starting 
from the fully aligned state |M rem (0, T, t = 0)| = 0.8, 
the remanent film magnetization relaxes fast to its equi- 
librium value. For the assumed temperatures the dipole 
coupling leads to a net magnetization |M rem | > |M r ° om | 
where for |M r ° cm | the dipole interaction is neglected. Af- 
ter several hundred MCS equilibration is obtained for the 
dipole-coupling induced |M rem | which then stays stable 
within the simulation time. We emphasize that it is im- 
possible to obtain these and the following results with 
conventional SSF algorithms. 

In Fig. I3b), the equilibrium magnetization 
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|M cq (0,T)| is shown as function of temperature T 
for different coverages 0. For low temperatures, clearly 
a magnetic ordering due to the dipole interaction is seen. 
The larger the coverage is, the larger is the ordering 
effect, since with an increasing the average island 
cluster size and thus the average dipole coupling energy 
increases. Above the ordering temperatures Tc(0), the 
magnetizations |M oq (0,T)| reach the corresponding 
values |Mg q (0,T)| as calculated without the dipole in- 
teraction. Due to the use of absolute values, |M° q (0, T)\ 
stays always finite. Tc(0) is estimated by extrapolating 
the linear part of |M eq | to |MgJ. The ordering temper- 
ature even for the largest investigated coverage is quite 
small, yielding T c « 6 K for = 0.8 ML. The rounding 
of |M oq | near the ordering temperature is caused by (i) 
the finite unit cell size, (ii) the use of the absolute value, 
(hi) the presence of island size and -position dispersions, 
and (iv) the average over 20 different realizations of the 
unit cell. 

The existence of a long-range magnetic ordering due 
to the dipole interaction has been calculated for peri- 
odic lattices^ In our study, we find that also within an 
irregular island system below the percolation threshold 
the dipole interaction leads to a magnetic ordering, in- 
dicated by a net magnetization |M oq | > |M° q |. Such a 
collective state is expected to be spin-glass-like, as dis- 
cussed for three-dimensional magnetic particle systems, 14 
and which needs a further investigation. We remark 
that the obtained ordering temperatures will increase by 
considering a noncollincar island magnetization beyond 
Sf = ±lj2& by taking into account the finite island ex- 
tension for the dipole interaction beyond the point-dipole 
approximation^ or by consideration of densely packed 
three-dimensional particles. 

In Fig. |31 we investigate the influence of the magnetic 
anisotropy and the exchange interaction on the magnetic 
relaxation for coverages < 0p. Here, the dipole cou- 
pling is neglected. In Fig. |3fa) , the remanent magnetiza- 
tion | M rem (0, T, t) | is given as function of MC time t for 
different temperatures T and for K = 0.01 meV/atom. 
The coverage is assumed to be = 0.8 ML. Starting 
from the fully aligned state, at first |M rem (0, T, t)\ drops 
rapidly due to relaxation of single islands and small is- 
land clusters. The further relaxation happens much more 
slowly since here larger island clusters have to be re- 
versed. For T > 25 K, the magnetization |M rcm (0, T, t)\ 
reaches within the depicted time range the curve |M r ° em | 
as calculated for K = 0. 

In Figs. Efl>) and Etc), we show |M rero (0, T, t)\ af- 
ter t — 1000 MCS for different coverages as function 
of temperature, using the anisotropy parameters K 
( ) .01 and 0.1 meV/atom. With increasing temperature, 
the magnetization \M mm (Q,T, t)\ approaches the equi- 
librium value |Afg q |. The corresponding blocking temper- 
atures Tb(0, K) are obtained by extrapolating the linear 
part of |M rem | to |M° |. A rounding of |M rem | near It, 
is observed due to the same reasons as discussed in con- 
nection with Fig. [21 



We emphasize that for a connected island structure an 
increase of the anisotropy K by a factor of 10 does not 
necessarily lead to an increase of T^(Q,K) by the same 
factor as obtained from the Stoner-Wohlfarth model. 19 
This is caused by internal cluster excitations, i.e., cre- 
ation or motion of domain walls inside island clusters. 
To discuss this, we have performed additional calcula- 
tions for an infinite domain wall energy 7, indicated by 
the full lines in Figs.|2tb) and|^c), hence allowing only 
for coherent island cluster rotations. Above a certain 
temperature, the curves for finite and infinite 7 deviate, 
since then internal cluster excitations become effective. 
For K — 0.01 meV/atom and coverage — 0.6 ML, the 
difference between these curves is small, thus the mag- 
netic relaxation happens mainly via coherent rotation. 
In contrast, for = 0.8 ML or for the larger anisotropy 
K = 0.1 meV/atom, obviously both relaxation processes 
are present. 

Which relaxation process is effective at a given temper- 
ature is determined by its energy barrier AE. For a co- 
herent rotation of an isolated island cluster, AE is given 
by its total anisotropy energy /C„. In contrast, AE for an 
internal cluster excitation consists of both the anisotropy 
of the actually reversed islands and the domain wall en- 
ergy, see Eqs.Q and JSJ. By closer investigation, we 
found that a particular relaxation process becomes ef- 
fective above a temperature amounting to 5 - 10 % of 
AE. In the temperature ranges T < 40 K for K = 
0.01 meV/atom and T < 100 K for K = 0.1 meV/atom, 
each internal cluster excitation consists mainly of revers- 
ing only one or two islands. For markedly larger temper- 
atures, the internal cluster excitations will become more 
complex, depending in a complicated way on the nanos- 
tructure. 

The influence of the dipole coupling on the relax- 
ation behavior is discussed in Fig. 0] For < 0p and 
K = 0.01 meV/atom, the dipole interaction results in a 
small increase of |Af rem | and It,. Thus, the dipole interac- 
tion leads to a slower magnetic relaxation. Interestingly, 
this effect is visible in the whole temperature range up 
to Tt>, and is not limited to those small temperatures 
where the dipole coupling induces a magnetic ordering, 
see Fig. [3] The increase of |M rem | will become larger if 
a stronger dipole coupling is assumed, for example for 
larger island magnetic moments. Experimentally, a sim- 
ilar effect was observed for two-dimensional arrays of in- 
teracting magnetic nanoparticles with random anisotropy 
axes^ A more detailed investigation of this property is 
needed. 

Next, we investigate the magnetization for coverages 
above the percolation threshold, > 0p. In Fig. [5] the 
equilibrium magnetization |M cq (0,T)| is shown as func- 
tion of temperature for different coverages. Here, the 
exchange coupling causes a fast magnetic relaxation and 
a strong ferromagnetic long-range order. For large cov- 
erages and at low temperatures, the behavior of |M eq | 
is reigned by the decrease of the internal island magne- 
tization TOi(0,T), whereas at elevated temperatures a 
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strong decay of |M eq | is caused by the disturbance of 
the island spin alignment. The resulting ordering tem- 
peratures To (6) are deduced from the inflection points 
of \M cq (Q,T)\. In addition, for 6 = 1.0 ML and 
1.2 ML we show the nonequilibrium remanent magne- 
tization |M r em| after t = 1000 MCS, considering K 
0.01 eV/atom. The corresponding blocking temperatures 
Tb(0) are markedly larger than Tc(0) in the coverage 
range > 0p, where the nanostructure of the perco- 
lated thin film is still very irregular and nonequilibrium 
effects due to anisotropy barriers are pronounced. For 
films with larger coverages, having a higher connectivity 
between islands, the exchange coupling results in a fast 
magnetic relaxation and thus the temperature difference 
between 7b and Tq is small. A very weak ordering effect 
due to the dipole interaction is visible only for very low 
temperatures and coverages > 0p. Here, still a few 
isolated islands or island clusters exist which are coupled 
to the percolating cluster by the dipole interaction. 

In Fig. |SJ we summarize the most important results 
of this study. The (nonequilibrium) blocking tempera- 
ture Tb(0,^O and the (equilibrium) ordering tempera- 
ture Tc(0) are presented as functions of coverage in 
the whole investigated growth range. For a better vi- 
sualization, a logarithmic temperature scale is applied. 
Tb is determined for two different anisotropies K = 0.01 
and 0.1 meV/atom and t = 1000 MCS. Below the perco- 
lation coverage 0p, the dipole interaction induces small 
ordering temperatures Tq of the order of 1 - 10 K for the 
assumed model parameters. Due to the coagulation of is- 
lands with increasing coverage the exchange interaction 
becomes more important, since it couples single islands 
to magnetically aligned large clusters. This results in a 
strong increase of Tq in particular close to 0p. This be- 
havior has been observed in experiments on Co/Cu(001) 
ultrathin films ( l Tc-jump")£& For percolated thin films, 
the ordering temperature is of the order of 100 - 300 K 
and is, within the accuracy of our calculations, exclu- 
sively determined by the exchange coupling. The slope 
of 2c(0) for > 0p is not as steep as for < 0p. In 
addition, we show the ordering temperature Tq due to 
dipole interaction by neglecting the exchange coupling 
between islands. Evidently, a distinct variation of Tq 
near 0p is not obtained in this case. 

The nonequilibrium behavior as caused by the 
anisotropy K differs strongly for coverages below and 
above 0p. Due to the slow relaxation of the irregular 
atomic structure for < 0p, a blocking temperature 
7b(0) is obtained which is an order of magnitude larger 
than Tc(0) resulting from the dipole interaction. Evi- 
dently, Tb depends on the anisotropy K and the MC time 
t. On the other hand, for > 0p, the relaxation is ac- 
celerated by the exchange interaction. With increasing 
the remanent magnetization reaches the equilibrium 
value within t = 1000 MCS, hence Tb(0) merges into 
T C (0). 

Recently, a mean field theory (MFT) for the dipole- 
coupling induced magnetic ordering temperature has 



been performed, using a simplified growth modelA A 
qualitatively similar behavior of Xc(0) as compared to 
the present MC calculations was obtained, in partic- 
ular the strong variation of Tq near 0p due to the 
exchange interaction between coagulated islands. Evi- 
dently, MFT yields much larger values for Tc(0) for such 
low-dimensional systems especially for < 0p, since 
thermal fluctuations are neglected. 

In the following, we discuss our results in connec- 
tion with measurements on Co/Cu(001) ultrathin films. 
Although several model parameters are chosen in ac- 
cordance with this system, a full quantitative compar- 
ison cannot be drawn yet. The main reason is that 
the observed intermixing of Co adatoms with Cu sub- 
strate atoms is not taken into account within our growth 
model due to the incomplete knowledge of the resulting 
atomic morphology. The measured percolation threshold 
0p w 1.7 ML is much larger than the one as obtained 
with the growth parameters used by us. 

To investigate solely the effect of an enlarged 0p, we 
have performed additional simulations simply by taking 
into account magnetic islands with up to three atomic 
layers, yielding the observed 0p. Then for a coverage 
= 1.6 ML, the dipole coupling induces a ordering tem- 
perature Tq ~ 50 K. The corresponding blocking tem- 
perature for K = 0.01 meV and t = 1000 MCS is ob- 
tained to be Tb ~ 150 K. Hence, we find that for the as- 
sumed growth modes and magnetic parameters the block- 
ing temperatures are always markedly larger than the 
ordering temperatures. These temperatures are compa- 
rable with measured temperatures Tc ~ 150 K for cover- 
ages slightly below 0p A In addition, to draw a quantita- 
tive comparison with the Co/Cu(001) system a fourfold 
symmetry of the in-plane anisotropy has to be taken into 
account. We expect that the above stated general be- 
havior of the magnetic relaxation and ordering obtained 
with a uniaxial anisotropy will not be changed. 



IV. CONCLUSION 

In this study, we have calculated the nonequilibrium 
and equilibrium remanent magnetization of growing ul- 
trathin films, using a cluster Monte Carlo method. An 
island-type nanostructure with a nonuniform distribu- 
tion of island sizes, shapes, and locations was investi- 
gated. Within a micromagnetic model (modified Ising 
model) the single-island magnetic anisotropy, the dipole 
coupling, and the exchange interaction between magnetic 
islands were taken into account. We have analysed the 
transition from dipole-coupled islands for film coverages 
below the percolation threshold 0p towards a connected 
ferromagnetic film above 0p with increasing film cover- 
age. 

For coverages < 0p, the dipole interaction leads 
to an equilibrium net magnetization refering to a collec- 
tively ordered state. A small ordering temperature Tq of 
about 1 - 10 K results for the assumed model parameters. 
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Tc increases strongly near Op due to exchange interac- 
tion which aligns coagulated islands. On the other hand, 
the anisotropy induces a pronounced nonequilibrium re- 
manent magnetization which may be visible in experi- 
ment even after long waiting times. The corresponding 
blocking temperature Tb is obtained to be of the order of 
10 - 100 K, which is always markedly larger than To Ap- 
proaching Op , the proportionality between blocking tem- 
perature and anisotropy is no longer valid due to relax- 
ation via internal island cluster excitations. A nonequi- 
librium remanent magnetization due to anisotropy is vis- 
ible also for coverages O > Op where the film is still 
very irregular. For smoother films at larger coverages the 
exchange interaction induces a fast magnetic relaxation 
towards equilibrium. 

We have obtained these results with a cluster-spin-flip 
algorithm which takes into account coherent magnetic ro- 
tations of island clusters. This method leads to a very 
fast and more realistic magnetic relaxation towards equi- 
librium in the coverage range with an irregular nanos- 
tructure. Our results cannot be achieved by conventional 
single-spin-flip algorithms. The suggested CSF algorithm 
can be applied also to other inhomogeneous spin systems 
such as diluted magnets and spin glasses. 



Several possible improvements of our micromagnetic 
model are pointed out. In this study, we have used Ising- 
like states Sf = ±1. By applying continuously vary- 
ing spins Si, noncollincar magnetic arrangements can be 
analysed^ allowing one to also determine the effects of 
an external magnetic field for these strongly inhomoge- 
neous films. In particular, the movement of magnetic 
domain walls can be investigated. Furthermore, various 
magnetic nanostructures like chains and stripes 3 - can be 
studied easily by a proper variation of the parameters of 
the Eden-type growth model. Anisotropies with, e.g., a 
four-fold in-planc symmetry will be considered. Finally, 
the relaxation laws and times of the remanent magneti- 
zation can be investigated for such thin film systems. 3 
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In addition to the magnetic dipole coupling other long- 
range interactions can be present such as the indirect ex- 
change (RKKY-) interaction in case of metallic substrates, 
and the superexchange for insulating substrates. 
Consider two disk-shaped Co islands with N atoms each. 
Assume that the distance R between their centers is com- 
parable to their diameters, R « D oc a yN. Then 
the point-dipole energy of this island pair is -Edip ~ 
{Nfi at ) 2 /R 3 » 0.623 {pt % /ai)VN ~ 6 K, with /x at = 
2.0 hb the atomic magnetic moment measured for a 2-ML 
Co film)— /is the Bohr magneton, a — 2.5 A the Co in- 
teratomic distance, and N = 1000. 

Consider for instance the magnetic rotation ff — H of 
an isolated island cluster consisting of two connected is- 
lands. Assuming Ising-like magnetic moments, inevitably 
the single-spin-flip method applies a subsequent rotation 
of the single island spins through the intermediate state 
|L However, such a process is very unlikely due to the 
large increase of exchange energy of this intermediate state. 



Obviously, a proper treatment of magnetic relaxation re- 
quires the inclusion of a coherent rotation of the island 
pair. This is performed within a cluster-spin-flip algorithm 
which takes into account such simultaneous spin rotations 
of connected islands. 

The described CSF algorithm satisfies the condition of de- 
tailed balance. This is guaranteed by the fact that the 
probabilities for construction and choice of the cluster C v 
are the same for both flip directions, and by the used flip 
rates which already obey detailed balance. Ergodicity is 
maintained since any spin state can be reached due to the 
allowance of single-spin flips. 

A constant attempt frequency r o is widely applied in liter- 
ature. The justification of this approach is given in Refs. 7 
and 8 where r o is calculated using different approaches. Y a 
is found to vary weakly with temperature and local fields. 
Note that the applied transition rates depend only weakly 
on the exact value of r o . 
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FIG. 1: Comparison of the calculated film magnetization using the single-spin-flip (SSF) and cluster-spin-flip (CSF) algorithm. 
Only the exchange coupling between magnetic islands is considered. The remanent magnetization |Af rom | as function of MC 
time t for two film coverages (a) = 1.8 ML and (b) O — 0.8 ML above and below the percolation coverage ©p ~ 0.9 ML is 
shown. Within the CSF method, different maximum numbers A max of coherently flipping island magnetic moments are used. 
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FIG. 2: Long-range magnetic ordering due to dipole interaction for coverages O below the percolation coverage ©p. Only 
dipole and exchange interactions are included. The quantities |MrJ, m j and |Me q | neglect the dipole interaction, (a) Remanent 
magnetization |M rom | as function of MC time t for = 0.8 ML < 0p and different temperatures T. (b) Equilibrium 
magnetization |M cq | as function of temperature T for different coverages < 0p. 
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FIG. 3: Relaxation of the remanent magnetization for coverages O below the percolation coverage Op. Only anisotropy and 
exchange interaction are taken into account, (a) Remanent magnetization |M rem | as function of MC time t for O = 0.8 ML < Op, 
anisotropy K — 0.01 meV/atom, and different temperatures T. The full line refers to K = and T = 35 K. Furthermore, 
(b) and (c) show the remanent magnetization |M rcm | after t = 1000 MCS as function of temperature T for K = 0.01 and 
0.1 meV/atom. The full lines are calculated with an infinite domain wall energy 7. 
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FIG. 4: Influence of magnetic anisotropy and dipole coupling on the remanent magnetization |M rom (T)| as shown in Fig.^b). 
Results for two coverages O below the percolation coverage Op are presented. 




FIG. 5: Long-range ferromagnetic ordering due to exchange interaction for different coverages above the percolation coverage 
Op. The equilibrium magnetization |M eq | is shown as function of temperature T. For comparison, for Q = 1.0 and 1.2 ML 
also the remanent magnetization |M rcm | after t = 1000 MCS is depicted for an anisotropy K = 0.01 meV/atom. 
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FIG. 6: Semi-logarithmic plot of the magnetic ordering temperature Tc and the blocking temperature Tb as functions of the 
film coverage O. The temperatures Tc and T b are extracted from the preceding figures. The entire investigated coverage range 
below and above the percolation coverage Op is shown. The different magnetic interactions dominate in different coverage and 
temperature ranges as indicated. The dashed curve refers to the dipole-coupling-induced ordering temperature neglecting the 
domain wall energy 7 between islands. 



